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The effect of morphologic factors on magnetic flux dynamics and critical currents in percolative 
superconducting structures is considered. The superconductor contains the fractal clusters of a 
normal phase, which act as pinning centers. The properties of these clusters are analyzed in the 
general case of gamma-distribution of their areas. The statistical characteristics of the normal 
phase clusters are studied, the critical current distribution is derived, and the dependencies of the 
main statistical parameters on the fractal dimension are found. The effect of fractal clusters of a 
normal phase on the electric field induced by the motion of the magnetic fiux after the vortices have 
been broken away from pinning centers is considered. The voltage-current characteristics of fractal 
superconducting structures in a resistive state for an arbitrary fractal dimension are obtained. It is 
found that the fractality of the boundaries of normal phase clusters intensifies magnetic fiux trapping 
and thereby increases the current-carrying capability of the superconductor. 

I. INTRODUCTION 

A significant property of clusters of a normal phase embedded in a superconducting medium consists in tlicir 
capability to trap a magnetic flux and hold the vortices from moving under the action of the Lorentz force. As a 
result, these clusters can act as effective pinning centers Ql- M. This property is widely used in creating composite 
superconducting materials with high critical current values Tsl , B . The morphologic characteristics of the normal phase 
clusters essentially affects the dynamics of trapped magnetic flux, especially if the clusters have fractal boundaries 0- 
1^ . In the present work we consider in detail the geometric probability properties of these fractal clusters and their 
influence on the critical current and dynamics of trapped magnetic flux near the transition of the superconductor into 
a resistive state. 

II. FRACTAL GEOMETRY OF NORMAL PHASE CLUSTERS AND DISTRIBUTION OF CRITICAL 

CURRENTS 

Let us consider a superconductor containing fragments of a normal phase. Suppose that dimension of these fragments 
along one direction signiflcantlv exceed the other sizes. Similar columnar defects are of great interest in creating 
artificial pinning centers [||, [|lO[- (l2). If such a superconducting structure is cooled in a magnetic field directed along 
the axis of alignment of these defects below the critical temperature, then the distribution of magnetic flux trapped 
in the clusters of a normal phase will be two-dimensional. This can be done especially easily with a superconducting 
film where such clusters are formed near defects at the boundary with the substrate during the growth process and 
are directed transversely to the film plane [|j , . Let us suppose that the fllm surface fraction covered by the normal 
phase is below the percolation threshold for the transfer of magnetic flux (i.e., less than 50% for 2D-percolation 
In this case the fraction of the superconducting phase exceeds the percolation threshold, so there is a superconducting 
percolation cluster that can carry a transport current in the film plane. This kind of structure provides for the effective 
pinning and thereby raises the critical current, for the magnetic fiux is trapped in isolated clusters of the normal phase 
and the vortices cannot leave them without crossing the surrounding superconducting space. As the current increases, 
the trapped flux will remain unchanged until the vortices begin to break away from those clusters for which the 
pinning force is less than the Lorentz force caused by the transport current. When the magnetic flux is breaking away 
from the pinning centers, each vortex must cross the infinite superconducting cluster. In this case the vortices will 
move primarily along the weak links that connect the clusters of a normal phase to each other ||], These 
weak links form readily in high-temperature superconductors (HTS) characterized by a small coherence length 
I p^ . Structural defects, which could act simply as scattering centers at a large coherence length, create weak links in 
HTS. There is a large variety of weak links in a wide range of spatial scales in HTS ||], IQ. At an atomic level 



1 



weak links are created by structural point defects, primaxily by oxygen vacancies |16|, il9(|. On a mesoscopic scale 
weak links are efficiently formed on twin boundaries , ||l6[ , , | 
tens of nanometers from each other, therefore even a single crystal can have a fine substructure created by twins. 
Finally, on a macroscopic scale weak links are created by various structural defects, such as boundaries of grains and 
crystallites or barriers arisin g fr om the secondary degrading of a non-stoichiometric crystal into the domains with a 
high and low oxygen content |T7|]- JTot . Moreover, a magnetic field favors the coherence length to further decrease plf , 
so weak link formation is facilitated still more. In conventional low-temperature superconductors, which are of large 
coherence length, weak links are formed due to the proximity effect at sites where the distance between neighboring 
clusters of a normal phase is minimum. 

Thus, regardless of their origin, weak links form channels for vortex transport where the pinning force is less than 
the Lorentz force caused by the transport current. According to their configuration each cluster of the normal phase 
has its own value of the depinning current that makes the contribution to the overall statistical distribution of the 
critical currents. When a transport current is passed through the sample, vortices will be broken away first from 
clusters of a smaller pinning force and, accordingly, a smaller critical current of depinning. So the change in the 
trapped magnetic flux will be proportional to the number of the normal phase clusters which have the critical current 
less than the preset value /. Therefore, the relative decrease in the flux can be expressed through the function of 
cumulative probability F = F (1) for the distribution of critical currents of the clusters: 

-^-FW (1) 

where 



20 1 . Twins can be located at a distance of several 
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The right-hand side of Eq. (|l|) is the probability that any j-th cluster has a the critical current Ij less than the 
given upper bound /. On the other hand, the magnetic flux trapped in a single cluster is proportional to its area 
A. Therefore, the decrease in the trapped flux can be expressed through the function of cumulative probability 
W = W (A) for the distribution of the areas of normal phase clusters, which is proportional to the number of clusters 
of area smaller than a given value of A 

-^ = l-WiA) (2) 

where 



W {A) = Pr {VAj < A} 

Generally, the distribution of cluster areas can be described by the gamma-distribution, for which the function of 
cumulative probability has the form 

WiA) = ^(s + '^^/^o) (3) 

where 7 {v, z) is the incomplete gamma function, T [v) is Euler gamma function, and Aq and g are the parameters of 
the gamma-distribution that determine the mean area of the normal phase cluster A = [g -\- V)Aq and its standard 
deviation a a — A^^^g + 1. The distribution function of cluster areas can be found from the geometric probability 
analysis of electron photomicrographs of superconducting films For example, in the practically important 

case of YBCO films with columnar defects the exponential distribution is realized, which is a special case of the 
gamma-distribution of Eq. with 5 = 0. 

To clarify how the transport current affects the trapped magnetic flux, it is necessary to find the relationship 
between the distribution of critical currents of clusters of Eq. (Q) and the distribution of their areas of Eq. §). The 
larger the normal phase cluster size, the more weak links rise from its perimeter bordering with the surrounding 
superconducting space, and therefore, the smaller is the critical current at which the magnetic flux breaks away 
from this cluster. Let us suppose that the weak link concentration per unit length of the perimeter is the same 
for all clusters and all the clusters of equal perimeters have the same pinning force. Then the critical current / is 
inversely proportional to the perimeter P of the normal phase cluster: / oc because the probability to find a 

weak link over the larger perimeter is higher. It is suggested here that the vortex, driven by the Lorentz force, passes 
through the weak link from one normal phase cluster to another with the probability of 100%. In this case magnetic 
flux is transferred through the superconducting medium by Josephson vortices. The probability to trap the vortex 
when it moves through the weak link under the action of the Lorentz force is negligible, inasmuch as the Josephson 
penetration depth in the considered materials significantly exceeds the sizes of possible irregularities along the weak 
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link. In accordance with these assumptions, in order to find the distribution function of Eq. we have to know the 
relationship between the perimeter and the area of a cluster. As was first established in Ref. ||^, fractal properties 
of the boundaries of normal phase clusters essentially affect the dynamics of magnetic fiux in superconductors. For 
fractal clusters the perimeter-area relation has the form 

P oc (4) 

where D is the fractal dimension of the cluster boundary. 

Relation of Eq. (^) agrees with the generalized Euclid theorem ||2^ , according to that the ratios of the corresponding 
measures are equal when they are reduced to the same dimension. So P^l^ cx A^/^, and this relation holds true both 
for Euclidean clusters (for which the Hausdorff-Bezikovich dimension of the perimeter is equal to the topological 
dimension of a line, D = 1), and for fractal clusters (for which the Hausdorff-Bezikovich dimension of the boundary 
strictly exceeds the topological one, £>>!). 

Note that it is just the statistical distribution of cluster areas, rather than their perimeters, is fundamental for 
deriving the critical current distribution. Inasmuch as the Hausdorff-Bezikovich dimension of a fractal line exceeds 
unity, the perimeter of a fractal cluster is not well defined: it diverges when the precision of its measurement increases 
indefinitely |^^. At the same time, the topological dimension of the cluster area coincides with the Hausdorff- 
Bezikovich dimension (both are equal to 2). Therefore, the area of a surface confined by the fractal curve is a finite, 
well-defined quantity. 

Analyzing the geometric characteristics of normal phase clusters, we are considering the cross-section of extended 
columnar defect by the plane the transport current flows through. Therefore, though a normal phase cluster is a self- 
affine fractal ||2^ , we can consider its geometric probability properties in the planar section only, where the boundary 
of cluster is statistically self-similar. 

Now, using relation of Eq. (^) between the fractal perimeter and area, together with the assumption of inverse 
proportionality of the critical current of a cluster to its perimeter, we arrive at the following expression for critical 
current of the cluster: / = aA^^/"^ , where a is a form factor. Taking into account our initial relations of Eqs. (|l|) and 
(H), we can find the distribution of critical currents in the general case of the gamma-distributed cluster areas: 

^ ^ r(5 + 1) ^ ^ 



where 



6(3+1 _ (15/2) exp (61) r(.g + 1,61); ' - • • ^ 

Tiy^z) is the complementary incomplete gamma function, i = I/Ic is the dimensionless electric current, and Ic = 

a (AqG) is the critical current of the transition into a resistive state. The found function of cumulative probability 
of Eq. (g) provides a complete description of the effect of the transport current on the trapped magnetic fiux. Using 
this function, we can easily derive the probability density f {i) ~ dF/di for the distribution of dcpinning currents, 
which has the form 

m = DrTg + if "^'^'^'' {-G^-i) (6) 

The probability density is normalized to unity over all possible positive values of the critical currents. The relative 
change in the trapped fiux A(f>/$, which can be directly found from Equation (||), also determines the density of 
vortices n broken away from pinning centers by the current of a given magnitude i 

B r B A* 

„,)._^/,.V.'.-- (7, 

where B is the magnetic field, and <i>o = he/ (2e) is the magnetic flux quantum [h is Planck constant, c is the velocity 
of light, and e is the electron charge). 

Figure |l| shows how g-parameter of the gamma-distribution affects the distribution of the critical currents of the 
clusters. As an example, we have taken a value of the fractal dimension D = 1.5, which is near to the value 
D = 1.44 ±0.02 found earlier |^ for the normal phase clusters in YBCO fllm structures. On the other hand, the value 
D — 1.5 is intermediate between two limiting cases: D = 1 for the Euclidean clusters and D = 2 for the clusters of the 
most fractality. When the fractal dimension differs from unity so much, the fractal properties of the cluster structure 
of a superconductor are of great importance. In Fig. ^ the cluster area distribution is presented. The corresponding 
probability density has the following form 
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^(")= rC9+Tr '^"^''''^° (8) 



where a = A/ A is the dimensionless cluster area, with a = 1 and a a = l/Vff + 1- The function of cumulative 
probability for the dimensionless area is related with function of Eq. (||) by the formula W (a) = w (a') da' and can 
be written as 

7((g + l,(ff+l)a) 

"^^"^^ rUTT) 

As is seen from Fig. |], when ^-parameter decreases, the distribution f — f if) spreads to the right, spanning 
higher and higher range of critical currents. It might be expected that the highest current-carrying capability of the 
superconductor can be achieved in the limiting case of g = 0, when the gamma-distribution of Eq. ^ is reduced 
to the exponential one w (a) — exp(— a). In this case the small clusters give the most contribution into the overall 
distribution (curve 1 in Fig. ^). For g — Q expressions of Eqs. (||) and (H) can be simplified 

F{e)^e^^(-m^Y^\-i\ (9) 



.--exp(^-(^^j r^j (10) 

where, as above, i = I /Ic, and the critical current of the resistive transition can be calculated using a simpler formula: 

h = {2/(2 + D)f+''^'^aA-°'\ 

Figure shows how the fractal dimension of the cluster boundaries affects the distribution of the critical currents of 
Eq. (IQ). ^s is clearly seen from the figure, with increasing fractal dimension the critical current distribution spreads 
and shifts to the higher values of current. This shift can be characterized by the dependence of the mean critical 
current on the fractal dimension, as is shown in Fig. ^. Although the mode of the distribution modef (i) — {2 + D) /2 
is linear with respect to the fractal dimension, the mathematical expectation obeys a much stronger superlinear law 
controlled by Euler gamma function 

2+ZJ 

' "''l-?) (11) 



2^ 

Thus, the geometric probability properties of the normal phase clusters determine the statistical distribution of 
depinning currents. Using this distribution along with Equation (|l|), the relative change in the trapped flux under 
the action of the transport current can be found. 



III. DYNAMICS OF MAGNETIC FLUX TRAPPED IN NORMAL PHASE CLUSTERS WITH FRACTAL 

BOUNDARIES 

The effect of the transport current on the trapped magnetic flux is demonstrated in Fig. ^. The change in the 
trapped flux was calculated using formula of Eq. (|l]) for the exponential-hyperbolic critical current distribution of 
Eq. (^. Curves 1 and 5 in Fig. [^correspond to the limiting cases of the Euclidean clusters {D = 1) and the clusters 
of maximum fractality (D = 2), respectively; they confine the region of changes in the trapped flux for all possible 
values of the fractal dimension. The function F — F [i) reveals one important property; namely, it is very flat near 
the coordinate origin. It can readily be shown that all its derivatives become zero at the origin point: d^F (0) /di^ = 
for any values of k. Therefore even the Taylor series expansion of the function in the vicinity of the origin converges 
to zero value but not to the quantity F itself. This mathematical feature has a clear physical sense: such a small 
transport current does not affect the trapped magnetic flux, because there are no pinning centers with such small 
critical currents in the total statistical distribution. The decrease in the magnetic flux becomes noticeable only near 
the point of a resistive transition (i = 1). A use of the exponential-hyperbolic distribution of critical currents of Eq. 
(9) excludes the uncertainty caused by truncation of non-physical negative values of the depinning currents, which 
may happen, e.g., in the case of the normal distribution [p5| , [2^ . 

As is seen from Fig. the breaking of the vortices away is observed mainly at i > 1, when the sample undergoes 
a transition into a resistive state. Figure demonstrates the practically important property of a superconducting 
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structure with fractal normal phase clusters: the fractality favors magnetic flux trapping, and thereby increases the 
critical current magnitude below which the sample remains in a superconducting state. Indeed, the transport current 
of the value i = 2 causes the 43% of the total trapped magnetic flux to break away from the usual Euclidean clusters 
{D = 1, curve 1), whereas this value equals only to 13.5% for the fractal normal phase clusters of the greatest possible 
fractal dimension {D = 2, curve 5). This is equivalent to the pinning increase of 218% in the latter case. The 
enhancement of pinning due to fractality can be characterized by the pinning gain factor 

A$(L» = 1) 

which is equal to the relative decrease in the fraction of the magnetic flux broken away from the fractal clusters of 
dimension D in comparison with the Euclidean ones {D = 1). In Fig. ^ the dependencies of the pinning gain factor 
on the transport current and on the fractal dimension are presented. The maximum pinning gain is achieved when 
the boundaries of clusters have the greatest possible fractal dimension {D = 2). Note that the pinning gain factor 
characterizes the properties of the superconductor in the range of transport currents corresponding to a resistive 
state {i > 1). With smaller currents, the trapped magnetic flux virtually does not change, because there are no 
pinning centers with such small critical currents (Figs. |l|, ^ and the breaking of the vortices away has not started 
yet. In the presence of a finite resistance any current flow is accompanied by energy dissipation. As for any hard 
superconductor (type-II, with pinning centers) the energy dissipation in a resistive state does not mean the phase 
coherence destruction yet. Some dissipation accompanies any motion of the magnetic flux; this effect can be observed 
in a hard superconductor even for small transport currents. Therefore, the critical current in such materials cannot 
be determined as the greatest dissipationless current. The superconducting state collapses only when the dissipation 
increases in an avalanche-like manner as a result of the thermo-magnetic instability. 

The reason for the pinning gain caused by the fractality of the cluster boundaries lies in the fundamental properties 
of the statistical distribution of critical currents (see Fig. In the case of Euclidean clusters the mean value of 
the critical current calculated in accordance with Eq. ( pT| ) is equal to i{D = 1) = (3/2)'^/^y^ = 3.2562, whereas for 
the clusters of the most fractality {D = 2) it is infinitely high. As is seen from Fig. ^, when the fractal dimension 
increases, the contribution of the clusters with high depinning currents to the overall statistical distribution increases, 
too, which results in the enhancement of magnetic flux trapping. 

In a resistive state a hard superconductor can be adequately described by its voltage-current characteristics. Using 
the fractal distribution of the critical currents of Eq. (||), we can find the electric field caused by magnetic flux 
motion after the vortices have been broken away from pinning centers. Inasmuch as each cluster of the normal phase 
contributes to the overall distribution of the critical currents, the voltage across the superconductor V = V (i) is a 
response to the sum of effects made by contribution from each cluster. This response can be expressed as a convolution 
integral 

i 

V = Rf j{i-i')J{i')di' (12) 



where i?/ is the flux flow resistance. This representation for the voltage across the sample is often used to consider the 
pinning of bunches of vortex lines in a superconductor |^8| , as well as to analyze the critical scaling of voltage-current 
characteristics |27|| , i.e., in all the cases when the distribution of depinning critical currents occurs. In the following 
analysis we will focus on the results of the properties of the exponential- hyperbolic distribution in Eq. (|l^); we will 
not consider questions associated with the possible dependence of the flux flow resistance Rf on the transport current. 

In the simplest case, when all pinning centers have the same critical current ic, all vortices are released simultaneously 
at i = ic, and their density, in accordance with Eq. (|^), has the form 

13 B 
n = — / 5 {i' — ic) di' — —h [i — ic) 
^0 Jo *o 

where S (i) is the Dirac delta function, and h{i) = ^^^^ i < ^^"^ Heaviside function. The trapped flux 

would change at once by 100% in this case: A<I>/<I> — h{i ~ ic). 

Thus, accordingly with Equation ([l^), for the (5- like distribution of the critical currents in the regime of viscous 
flow of the magnetic flux the voltage across the superconductor is controlled by a simple linear dependence V = 
Rf {i — ic) h(i — ic)- The corresponding voltage-current characteristic is shown in Fig. [t] by the dashed line a. In the 
same figure the voltage-current characteristics of the superconductor in a still simpler approximation of the critical 
state model is represented, too (the dashed curve b). In this case the response on any external action, which results 
in the appearance of an electric field in a hard superconductor, is the flow of the current that is equal to the critical 
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one i = ic, regardless of the value of the voltage across the sample. (The dimensionless critical current is equal to 
unity due to the normalization chosen above, i = I /Ic) 

For the fractal distribution of the critical currents the situation changes radically, because now the vortices are 
being broken away in a wide range of transport currents. We shall consider the case of the exponential distribution of 
the cluster areas {g = 0), where the current-carrying capability of the superconductor is maximum. After substitution 
of the distribution of Eq. ( |lO| ) into Eq. (^2|) followed by integration, the voltage across the sample can be expressed 
through the function of cumulative probability of Eq. (M) 



V Rf F{i')di 
Jo 



Upon integration we get 



V = Rf 



I exp 




(13) 



In the limiting cases oi D — 1 and D = 2, Equation can be simplified. 

For Euclidean clusters {D = 1) the voltage across the superconductor has the form 



V = Rf 



I exp 



3.375 



V3.3757rerfc 



(14) 



where the complementary incomplete gamma function is expressed through the complementary error function 
r(l/2,z) = V^erfc(Vi). 

For clusters of maximum fractality {D — 2), after substitution of the representation for the complementary incom- 
plete gamma function F (0, z) = — Ei (— z). Equation ( |l3| ) for the voltage takes its final form 



V = Rf 



I exp 



4Ei 



(15) 



where Ei(— z) is the exponential integral function. 

In Fig. 1^ there are shown the voltage-current characteristics of a superconductor containing fractal clusters of a 
normal phase. For all values of the fractal dimension there is a noticeable voltage drop starting from the transport 
current of i = 1 that coincides with the value found earlier from the critical current distribution of Eqs. (|^) and (|^) 
for the resistive transition current. Equations ( |l^ and (15) describe the dependencies of the voltage on the transport 
current for the limiting values of the fractal dimension. Whatever the geometric probability properties of the normal 
phase clusters may be, the voltage-current characteristics of a superconductor with these clusters will lie in the region 
confined by these dependencies (curves 1 and 5 in Fig. 0). As is seen from the figure, the fractality significantly 
reduces the electric field induced by magnetic flux motion. In Fig. || the dependencies of the attenuation factor of 
dissipation 
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V{D) 



dB 



on the transport current at various values of the fractal dimension are presented. The decrease in the electric field 
with increasing fractal dimension is especially appreciable in the range of currents 1 < i < 3, where the pinning gain 
also has a maximum (Fig. ph. Both of these effects have the same nature, since their reason lies in the peculiarities of 
the fractal distribution of the depinning currents. As is seen from Figs. H and ^ an increase in the fractal dimension 
leads to considerable spreading of the tail of the distribution f — f {i). This means that more and more clusters 
of small sizes, which can best trap the magnetic flux, are being involved in the game. As a result, the density of 
vortices broken away from pinning centers by the Lorentz force decreases, so the smaller part of magnetic flux can 
flow creating a lower electric fleld. In turn, the smaller the electric field, the smaller energy is dissipated when the 
transport current passes through the sample. So the decrease in heat evolution, which could cause a transition into 
a normal state, leads to an increase in the current-carrying capability of the superconductor containing such fractal 
clusters. 

Thus, Figs. clearly demonstrate the most important result: the fractality of the normal phase clusters prevents 
destruction of the superconductivity by the transport current, and thereby enhances the current-carrying capability 
of the superconductor. The fractal clusters essentially aflFect the dynamics of the magnetic flux trapped in the super- 
conductor. The crucial change of the depinning current distribution caused by increasing the fractal dimension of the 
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clustes is the reason of this effect. The fractahty of the boundaries of the normal phase clusters intensifies pinning and 
slows down destruction of the superconductivity by the transport current. This point provides principally new possi- 
bilities for increasing the critical currents of composite superconductors by optimizing their geometric morphological 
properties. 

This work is partially supported by Russian Foundation of Fundamental Researchs (Grant No 02-02-17667). 



[1] T. Higuchi, S. I. Yoo, and M. Murakami, Phys. Rev. B 59, 1514 (1999). 

[2] L. Krusin-Elbaum, G. Blatter, J. R. Thompson, D. K. Petrov, R. Wheeler, J. UUmann, and C. W. Chu, Phys. Rev. Lett. 

81, 3948 (1998). 

[3] S. N. Dorogovtsev and Yu. I. Kuzmin, Phys. Lett. A 170, 245 (1992). 

[4] R. Surdcanu, R. J. Wijngaarden, B. Dam, J. Rector, R. Griessen, C. Rossel, Z. F. Ren, and J. H. Wang, Phys. Rev. B 58, 
12467 (1998). 

[5] E. Mezzetti, R. Gerbaldo, G. Ghigo, L. Gozzelino, B. Minetti, C. Camerlingo, A. Monaco, G. Cuttone, and A. Rovelli, 
Phys. Rev. B 60, 7623 (1999). 

[6] M. R. Beasley, Artificially-Structured Superconductors. - In: Percolation, Localization and Superconductivity, Ed. by 
A. M. Goldman and S. A. Wolf, NATO AST Scries, Scr.B, vol.109, (Plenum Press, New York, 1984), pp.115-143. 

[7] C. J. Olson, C. Rcichhardt, and F. Nori, Phys. Rev. Lett. 80, 2197 (1998). 

[8] Yu. L Kuzmin, Phys. Lett. A 267, 66 (2000). 

[9] Yu. L Kuzmin, Tech. Phys. Lett. 26, 791 (2000). 
[10] L. Krusin-Elbaum, L. Civale, G. Blatter, A. D. Marwick, F. Holtzberg, and C. Feild, Phys. Rev. Lett. 72, 1914 (1994). 
[11] Yu. 1. Kuzmin and I. V. Plcshakov, Tech. Phys. Lett. 25, 475 (1999). 
[12] Yu. I. Kuzmin, I. V. Plcshakov, and S. V. Razumov, Phys. Solid State 41, 1594 (1999). 
[13] D. Stauffer, Phys. Rep. 54, 2 (1979). 

[14] V. V. Bryksin, A. V. Goltsev, S. N. Dorogovtsev, Yu. I. Kuzmin, and A. N. Samukhin, J. Phys.: Condens. Matter 4, 1791 
(1992). 

[15] M. J. M. E. de Nivelle, G. J. Gerritsma, and H. Rogalla, Phys. Rev. Lett. 70, 1525 (1993). 

[16] D. J. Scalapino, Phys. Rep. 250, 329 (1995). 

[17] R. Haslinger and R. Joynt, Phys. Rev. B 61, 4206 (2000). 

[18] H. R. Kerchner, D. P. Norton, A. Goyal, J. D. Budai, D. K. Christen, D. M. Kroegeer, M. Paranthaman, D. F. Lee, F. A. 

List, R. Feenstra, and E. H. Brandt, Phys. Rev. B 60, 6878 (1999). 
[19] A. I. Rykov, S. Tajima, F. V. Kusmartsev, et al., Phys. Rev. B 60, 7601 (1999). 
[20] I. Maggio-April, C. Renner, A. Erb, et al.. Nature 390, 487 (1997). 

[21] J. E. Sonier, R. F. Kiefi, J. H. Brewer, D. A. Bonn, S. R. Dunsiger, W. N. Hardy, R. Liang, R. I. Miller, D. R. Noakes, 

and C. E. Stronach, Phys. Rev. B 59, R729 (1999). 
[22] B. B. Mandelbrot, Fractals: Form, Chance, and Dimension (Freeman, San Francisco, 1977). 
[23] B. B. Mandelbrot, The Fractal Geometry of Nature (Freeman, San Francisco, 1982). 

[24] B. B. Mandelbrot, Self-affine Fractal Sets. - In: Fractals in Physics, Ed. by L. Pietronero and E. Tosatti (North Holland, 
Amsterdam, 1986), pp.3-28. 

[25] Yu. I. Kuzmin, A. P. Paugurt, I. V. Plcshakov, and S. V. Rasumov, Supercond. Sci. Technol. 7, 41 (1994). 
[26] R. Wordonwober, Phys. Rev. B 46, 3076 (1992). 
[27] B. Brown, Phys. Rev. B 61, 3267 (2000). 

[28] W. H. Warnes and D. C. Larbalestier, Appl. Phys. Lett. 48, 1403 (1986). 



7 



1.6 




FIG. 1. Distribution of the critical currents for various values of (;-parameter of the gamma-distribution at the fractal 
dimension of D = 1.5. Curve (1) corresponds to the case of g = 0, (2) - g = 1, (3) - g — 2, (4) - g — 5, and (5) - g = 10. 
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FIG. 2. Gamma-distribution of the areas of normal phase clusters. Curve (1) corresponds to y = 0, (2) - y = 1, (3) - y = 2, 
(4) - = 5, and (5) - = 10. 
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FIG. 3. Effect of tfie fractal dimension of tlie cluster boundary on the critical current distribution a,t g — 0. Curve (1) 
corresponds to the case of Euclidean clusters {D = 1), curve (2) - to the clusters of fractal dimension D = 1.25, {3) - D — 1.5, 
{4) - D — 1.75, and (5) - to the clusters of boundaries with maximum fractality {D — 2). 
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FIG. 4. Dependence of the mean critical current on the fractal dimension. 
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FIG. 5. Effect of a transport current on the magnetic flux trapped in fractal clusters of a normal phase at <; = 0. Curve (1) 
corresponds to D = 1, {2) - D = 1.25, (3) - D = 1.5, (4) - D = 1.75, (5) - £> = 2. 
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